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ABSTRACT 

In this paper, we study a Sturm-Liouville operator with eigenparameter- 
dependent boundary conditions and transmission conditions at three interior 
points. We give an operator-theoretical formulation, construct fundamental 
solutions and investigate some properties of the eigenvalues and corresponding 
eigenfunctions of the discontinuous Sturm-Liouville problem. 



1 Introduction 

In recent years, more and more researchers are interested in the discontinuous 
Sturm-Liouville problem for its application in physics (see [1-15]). Such prob- 
lems are connected with discontinuous material properties, such as heat and 
mass transfer, varied assortment of physical transfer problems, vibrating string 
problems when the string loaded additionally with point masses, and diffraction 
problems [8,9]. Moreover there has been a growing interest in Sturm-Liouville 
problems with eigenparameter-dependent boundary conditions, i.e., the eigen- 
parameter appears not only in the differential equations but also in the boundary 
conditions of the problems (see [1, 4-7, 10-13] and corresponding bibliography). 

In this study, we examine eigenvalues and eigenfunctions of the differential 
equation 

Tu :— ~a{x)u" (x) + q{x)u{x) = \u{x) (1) 
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on [a, hi) U {hi, /i2) U (/12, /13) U (/13, b] , with boundary conditions 

Tiu := aiM (a) + a2u'{a) = 0, (2) 
T2U := A + A) u (6) - (/3^A + P2) u'{b) = (3) 

and transmission conditions at the points of discontinuity x = hi (i = 1, 2, 3), 



T3U 


■= 71 (hi 


-0) 


- 6iu {hi + 0) 


= 0, 


(4) 




■■= I'l (hi 


-0) 


- S[u {hi + 0) 


= 0, 


(5) 


T5U 


■= 72 ih2 


-0) 


- 52U {h2 + 0) 


= 0, 


(6) 


Teu 


■■= 72 {h2 


-0) 


- {h2 + 0) 


= 0, 


(7) 




■■= 73 {ha 


-0) 


- 5su {hi + 0) 


= 0, 


(8) 




:=73(/i3 


-0) 


- S'^U {h3 + 0) 


= 0, 


(9) 



where a{x) = af, for x e [a, hi), a{x) = a^, for x G {hi,h2) ,a{x) = a| for 
X e {h2,h3), a{x) = a1 for x <E {hs,b], a < hi < h2 < h^ < b, q{x) is 
a given real- valued function continuous in [a, hi) , {hi, ^12) , (ft.2, ft-a) , (ft-a, and 
has finite limits q{hi ± 0) = \imx^h,,±o q{x) {i = 1, 2, 3); A is a complex eigen- 
value parameter; 7^, 5i,^[, 5[ {i = 1, 2, 3) , aj, a'j, 0j, f3j {j = 1, 2) arc real num- 
bers; |ai| + \a2\ + 0, |7i| + %\ and + %\ / = 1,2,3) and 

p = P'ip2 - PlP'2 > 0. 



2 Operator formulation and some basic proper- 
ties of this operator 

In the Hilbert space H = L2 {a, 6) © C we define an inner product by 



h2 hs 



{F, G):=^Jf {X) g{x)dx+^Jf {x) g {x)dx + ^^^^ J f i-) 9 {x)dx+ 



hi /i2 

3 

■ Si 5', 



for 



F := I "f^"! ) ■ G := ( S^^' 1 £ H. 



Following [7] for convenience wc put 

Ri {u) :=Piu{l)-P2u'{l), 
R[{u) := I3[u{l)-/3W{1)- 
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For functionals / (x) , which defined on [a, fti)U(/ii, /i2)U(/i2, ^3)U(/i3, b] and has 
finite hmits / {hi ± 0) := hm^^?,,±o / (x) {i = 1, 2, 3). By fi {x) {i = 1, 2, 3, 4) 
we denote the functions 

fi {x) := 
/2 (x) := 



,/3(x) := 



U ix) := 

which are defined on fli := [a, hi] , O2 := [/ii,/i2] ,^^3 := [h2,hs\ , O4 := [/is, 6] 
respectively. We can rewrite the considered problem (1) — (9) in the operator 
formulation as 

AF = XF 



fix), 




X G [-l,hi), 




^hi-ofix), 


X = hi, 


lima;_ 


^hi+ofix), 


X = hi, 




fix). 


X e (/ii,/i2) , 


lima;_ 


-^h2-ofix), 


X = /l2, 


lim^;^ 


-^h2+0 fix), 


X = h2, 




fix), 


X G (/l2,/l3) , 


limj:_ 


-^ha-afix), 


X = h3. 


lima;- 


^hs+ofix), 


X = hs. 


fix), 




X G {h3,b] 



where 



and with 



^^=( -Tif) '^""^^^ 



D (A) := i^F e H fiix) , f- (x) are absolutely continuous in fli (i = 1, 4) , r/ G 
L2 [a, b] , nf = T2f = ... = rs/ = and fi = R[ (/)} 



and 

rf 



AF 



-Ri if) 



Consequently, the problem (1) — (9) can be considered as the eigenvalue problem 
for the operator A. 

Theorem 2.1. The linear operator A is symmetric in H. 
Proof. Let F, G £ D (A). By two partial integrations, we get 

{AF, G) = {F, AG) + (W (/, g;hi-0)-W (/, g; a)) 

+ ^{W if, g;h2-0)-W if, g; hi+0)) 

in I 

+ iW if, 9; hs-0)-W {f,g; h^ + 0)) + 

7i7i7272 
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+ 



{R[{f)R,{g)-Rr {f)R'Ag)) 



(10) 



where 

W{f,g;x) = f{x)g'{x)-f'{x)gix) 

denotes the Wronskians of the functions / and g. Since / and g satisfy the 
boundary condition (2), it follows that 

W{f,g;a)=0. (11) 

Prom the transmission conditions (4)-(9) we get 

W if, g;h^-0) = ^W if, g; hi+Q), 



71 7i 

W if, g; h2-0) = (/, g; h2 + 0), 

3 

W{f,g;hs-0) = l[^W{f, g; + 0) . 



(12) 



i=l 



Furthermore, 

R[if)Ri {g)-Ri if)R'Ag) 

= i0[f (b) - P'2f' m m{h) - p^g' m - m (b) - m' (^>)) iP'igih) - nksib)) 

= (/32/3i - /3^/3i) /' (b) g (b) + - f (6) g' (b) 

= P if ib) gib)-f ib) g' ib)) - -pW if, g; b) . (13) 
Finally, substituting (11)-(13) in (10) then we get 

{AF,G) = {F,AG) iF,G&DiA)) (14) 

■ 

Corollary 2.1. All eigenvalues of the problem (l)-(9) are real. 

Wc' c;an now assume that all eigenfunctions are real- valued. 

Corollary 2.2. // Ai and A2 are two different eigenvalues of the problem 
(l)-(9), then the corresponding eigenfunctions ui and U2 of this problem satisfy 
the following equality: 

^ /ii 5 5'^^ 5 5' 5 5' ''^ 

ui (x) U2 ix) dx + „^ ^ , / ui (x) U2 (x) dx + „^ ^ ,^ ^ , / ui (x) U2 (x) dx 
aV ahil[J 0111111212 J 



1 1 lili 



/ ui (x) U2 (x) dx = -'-^ — R'l (Ul) R'l iU2) 
nt-,liliJ aip 



In fact this formula means the orthogonality of eigenf unctions ui and U2 in the 
Hubert space H. 

We need the following lemma, which can be proved by the same technique 
as in Theorem 1.5 in [15]. 

Lemma 2.1. Let the real-valued function q{x) be continuous in [a,b] and 
f {X) ,g (A) are given entire functions. Then for any A € C the equation 

—u" + q (x) u = Au, X G [a, b] 

has a unique solution u — u{x, A) satisfies the initial conditions 

u{a) = f (A) , u' (a) = g (A) (or u{b) = f (A) , u' (6) =g{\)). 

For each fixed x G [a, b] , u (x, A) is an entire function of A. 
We shall define two solutions 



01A {x), X G [a, hi) , 
S,(r\-l (t>2\{x), xeihi,h2), 
'^"^''^-^ hxix), xeih,,hs), ^^d^^W 

04A (x) , X e (/l3,6] , 



Xix {x), xe [a,hi) , 
X2\ (x), x G {hi,h2) ■ 

X3A (x), X € {h2, hs) , 
X4A (x), X G (/l3,6] , 



of the Equation (1) as follows: Let <j)ix {x) := <j)i {x, A) be the solution of Equa- 
tion (1) on [a, hi], which satisfies the initial conditions 

u (a) = a2, u' (a) = —ai. (15) 

By virtue of Lemma 2.1, after defining this solution, we may define the solution 
(j)2 {x, A) := (/>2A {x) of Equation (1) on [hi, /12] by means of the solution (j)i (x. A) 
by the initial conditions 

u {hi) = ]i</.i {hi. A) , u' {hi) = (/ii. A) . (16) 

Ol Ol 

After defining this solution, we may define the solution 03 {x, A) := (p^x {x) 
of Equation (1) on [/i2,/i3] by means of the solution 02 (2;, A) by the initial 
conditions ^ 

u {h2) = ]^02 {h2, A) , u' (/12) = {h2, A) . (17) 
62 b'2 

Continuing in this manner, we may define the solution 04 {x. A) := 04a {x) of 
Equation (1) on [/13, 6] by means of the solution 0^ {x, A) by the initial conditions 

u (/13) = ^03 {hs, A) , u' {hm) = ^(P's (/i3, A) . (18) 

Therefore, (a;, A) satisfies the Equation (1) on [a, /ii)U(/ii, /i2)U(ft-2, /i3)U(/i3, b], 
the boundary condition (3), and the transmission conditions (4)- (9). 

Analogically, first we define the solution X4A {x) Xi {x, A) on [/13, b] by the 
initial conditions 

u {b) = /3^A + /32, u' (6) = p'lX + (3i. (19) 
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Again, after defining this solution, we may define the solution X3X {x) ■= Xs {^^ ^) 
of the Equation (1) on [h2,h3] by the initial conditions 

u (/13) = -X4 (/is, A) , u' {hs) = 4x4 (hs, A) . (20) 

73 73 

Again, after defining this solution, we may define the solution X2X {x) '■= X2 {x, A) 
of the Equation (1) on [hi, /12] by the initial conditions 

u (/12) = -X3 {h2, A) , u' {h2) = 4x3 {h2, A) . (21) 
72 72 

Continuing in this manner, we may define the solution XiA {x) ■= Xi {x, A) of 
the Equation (1) on [a, hi] by the initial conditions 

u {hi) = ^X2 {hi, A) , u' {hi) = 4x2 {hi,X) . (22) 
71 7i 

Therefore, x {x, A) satisfies the Equation (1) on [a, /ii)U(/ii, /i2)U(/i2, /i3)U(/i3, b], 
the boundary condition (3), and the transmission conditions (4)-(9). 

It is obvious that the Wronskians 

uji (A) := Wx {(t>i,Xi\ x) := {x, A) x!i {x, X)-(l>'i {x, A) Xi {x, A) , a; G fi, (i = T74) 

are independent of x EVli and entire functions. 

Lemma 2.2. For each A G C, wi (A) = ^W2 (A) = :^^|^W3 (A) = 

Proof. By the means of (16)-(18) and (20)-(22), the short calculation gives 

Wx {^i,xi;hi) = (<^2, X2; hi) = ^Wx (.^3, X3; h2) = ^^^j^Wx {h, Xs; h2) 

7i7i 7i7i 7i7i7272 

= M^Wx{^„X4;hs)=(f[^]wx{<^„xr,hs) 

717W272 yfJi'yi'yiJ 

so wi (A) = 5li02 (A) = 5ld^oj3 (A) = dfS^dluji (A) . ■ 
Now we may introduce the characteristic function 

Lj (A) := uji (A) = 6I0J2 (A) = SjSluJs (A) = SfSldlcOi (A) . 

Theorem 2.2. The eigenvalues of the problem (l)-(9) are the zeros of the 
function lu (A) . 

Proof. Let uj {Xq) — 0. Then Wxq {(j)i,xi',x) = and therefore the functions 
<piXo {x) and x^Xo {x) are linearly dependent, i.e. 

XiAo {x) = ki(j)iXo {x), X G [a, hi] 
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for some ki ^ 0. From this, it follows that x (x, Ag) satisfies also the first 
boundary condition (2), so x (a;, Aq) is an eigenfunction of the problem (1) — (9) 
corresponding to this eigenvalue Aq. 

Now wc let ?io (x) be any eigenfunction corresponding to eigenvalue Aq, but 
u;(Ao) 7^ 0. Then the functions (/>i,Xi) '/*2, X2) '/>3) X3) '/*4) X4 would be linearly 
independent on [a, /ii] , [/ii, /12] , [/i2, /^s] and [/13, 6] respectively. Therefore uq (x) 
may be represented in the following form 



uq (x) 



ci(/)i (x, Aq) +C2X1 (a;, Aq) , xe[a, /ii), 

c3(/)2 (x, Ao) + C4X2 (a;, Ao) , a; e {hi,h2) , 

C5^3 (a;, Ao) + ceXs (a;, Aq) , a; e (/12, /13) , 

C7^4 {x, Ao) + C8X4 (a;, Ao) , x e (/13, 6] . 



where at least one of the constants ci, C2,...,cs is not zero. Considering the 
equations 

T„(uo(a;)) =0, u = T78 (23) 

as the homogenous system of linear equations of the variables ci, C2, cg and 
taking (16)-(18) and (20)-(22) into account, it follows that the determinant of 
this system is different from zero. Therefore, the system (23) has only the trivial 
solution a = (z = 1,8). Thus we get a contradiction, which completes the 
proof. ■ 

Lemma 2.3. // A = Aq is an eigenvalue, then (j){x,Xo) and x(a;, Aq) are 
linearly dependent. 

Proof. Let A = Aq be an eigenvalue. Then by virtue of Theorem 2.2 
W {(t)iXa , XiXo ;x)=u}i (Ao) = 

and hence 

X^Ao (x) = h(f),Xo {x) {i = M) (24) 

for some fci 7^ 0, fc2 ^ 0, ^ 0,k4 ^ 0, We must show that k\ = k2 = kz = k^. 
Suppose, if possible, that kz 7^ k^. 

Taking into account the definitions of the solutions 0, {x, A) and Xi {x, A) 
from the equalities (24), we have 

ri (XAo) = XAo (^^3 - 0) - — XAo {hz + 0) = X3Ao {hi) - — X4Ao {hi) 

73 73 

^3 *^3 ^3 

= kz(t>?, (/13) A;404 (/13) = ^3 — 04 (/13) A:404 (/13) 

73 73 73 

= — (fc3 - ki) <j>i {hs) = 0. 

73 

since T7 (xAq) = and ^ (^3 — k^) ^ 0, it follows that 

04Ao (hs) = 0. (25) 
By the same procedure from ts (xaq ) = we can derive that 

04Ao (hs) = 0. (26) 
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From the fact that (f)4x„ (x) is a sohition of the differential Equation (1) on [/13, b] 
and satisfies the initial conditions (25) and (26), it follows that (p^x^ (x) = iden- 
tically on [/13, b] because of the well-known existence and uniqueness theorem for 
the initial value problems of the ordinary linear differential equations. Making 
use of (16)-(18) and (25)-(26), we may also derive that 

hxo (hm) = (t)'zx„ (hs) = 
Continuing in this matter, we may also find that 

hxo ih2) = c^'sxo {h2) = 0. 
<p2Xo (hi) = ct>'^x, (hi) = 

<AiAo (^1) = ^1x0 ('^i) = 0- 

identically on [/12, ^3] , [fti, ft.2] , [a, /ii] respectively. Hence (/)(x, Ao) = identi- 
cally on [a, hi) U {h\,h2) U (/i2,/i3) U (/i3,6]. But this contradicts with (15). 
Hence ^3 = ki. Analogically we can prove that = fcs and ki = k2- ■ 

Corollary 2.3. If X = Xq is an eigenvalue, then both (j) (x, Aq) and x {x, Aq) 
are eigenf unctions corresponding to this eigenvalue. 

Lemma 2.4. All eigenvalues An are simple zeros of w(A). 
Proof. Using the well-known Lagrange's formula, it can be shown that 

(hi /l2 ft,3 

— j 4>X {x) (/>A„ {x) dx+ — (j)x (x) 0A„ (x) dx+— (l)x (x) 0A„ (x) dx 
a hi h2 



+ ^j<t>x {x) ct>x^ {x)dx \ =W{<f>x, ct>x„ ■,b) (27) 

for any A. Recall that 

Xx„ {x) = kn(f)Xn {x), X e [a, hi) U (/ii, /12) U (/12, /13) U (/13, 6] 

for some fc„ ^ 0, n = 1,2, .... Using this equality for the right side of (27), we 
have 

W{^x,<l^x^;b) = ^W{<l>x,Xx^;b) = ^{XnR[{<l>x) + Ri iM) 
= }-[u:{X)-{X-Xn)R'M] 



(A - A J ^ 



Ri [<px) 



X — Xn 

Substituting this formula in (27) and letting A — >^ A„, we get 

hi h2 /13 



(28) 



4 / (?^A„ (a;))^ dx + ^ [ ((/>A„ (x))^ <ix + ^ [ 



{<t>Xn {x))^ dx+\ I ((/)A„ (a;))^ dx 



2 

■4 

hi h2 /13 



= ^{^'{K)-K{ct>xJ). 
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Now putting 
in (28) we get w' (A„) ^ 0. 



3 Asymptotic approximate formulas of cu (A) for 
four distinct cases 

We start by proving some lemmas. 

Lemma 3.1. Let <p [x, A) be the solutions of Equation (1) defined in Section 
2, and let A = .Then the following integral equations hold for k — Q,l : 

Ah) , , / s {x — a) \ ai I . s [x — a) ^ 



4'i\ {x) = I COS I — ai — I sin 



. s{x-yy 



+ I ("sin " 5 (y) {y) dy, (29) 

aisj V ai J 

a 

W = ir^iA ih,) (cos ) + t ^0;, (h,) (sin £(2^) 

, _ (30) 

^^7^ j Q (y) 02A (2/) t^y, 



a2S 
h 



1 / / ■ SI 

— ' sm 



(fe) (31) 

9 (y) ?!'2A (y) ^2/, 



035 



-/•Ia (^) = if ^3A (/.3) (cos ^) + (/,3) (si 



Sill - 

a4 



+^/(«-^)^'^(2/)'^4A(y)dy, ^^'^ 
«;ftere(.)('=)=^(.). 

Proof. It is enough to substitute s^(f>i\ {y)+(j)'ix iv) ' ■*^'/'2A (y)+02A (2^) ' «^'/'3A (y)+ 
</'3A (y) ' «^'^4A (y) + iplx iv) instead of (y) (/)ia (y) , q {y) <l>2\ {y) , q (y) <l>3\ {y) , 
q (y) (j)4x (y) in the integral terms of the (29)- (32), respectively, and integrate by 
parts twice. ■ 

Lemma 3.2. Let X = s^, Ims = t. Then the functions dj\ (.7;) ha,ve the 
following asymptotic formulas for |A| — >■ 00, which hold uniformly for x & ili 
[for i = T~4 and k = 0, 1.) : 
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If a2j^0 



,w f s{x-a)Y''^ f k-2 |i|(a;-a) 
(t)ix [x) = a2 cos ^ +0[\s\ exp ' " 



(fe) 71/ s{x-hi)\^^ s{hi-a) 

92X [x) = a2-r cos cos 

Oi V 0-2 J ai 

a2li f . s (x - hi)\^'^'' . sihi-a) 
— 0L2 — — sm ■ 



a\b'i \ 0.2 ) 0'\ 

( fc_i |i| [ai (a; - /ii) + a2 (/ii - a)] 

V 01^2 



, , 7172/ s{x-h2)\^^^ s{h2-hi) s{hi-a) 

^3A (x) = 0C2^^ ( COS ) COS ■ ' 

(fe) 



5i52 V as 
7i72 / s{x - /12) 



0-2 -Frr- COS sm sm 

5[52 V «3 ' 

717^ ( . s{x-h2) Y''^ , _ 
a2 -r-rj- sm sm cos 

dido V "3 



7W2 f- s{x- /12) 



3 

a2 




S (/l2 - 


hi) 


02 




s {h2 - 


hi) ^ 


a2 




S (/l2 - 


hi) 



3 

Oi 




s {hi - 


a) 


ai 




.S {111 - 


a) 


ai 




s{hi - 


a) 



(fe) 

- Ct2 TTTT 1 sm ^ — I COS — sm 

O1O2 V «3 / (12 ai 

+ 0(|s|*~^ex 1^1 ["1^2 {x - h2) + aitts {h2 - fei) + a2a3 {hi - a)] 

010203 
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(k) 717273 ' O.^—M^InW 



(Pa\ yx) = a2 f c c COS COS cos 

71727302 / s{x - h:i) \ s (/13 - /i2) . s (/i2 - hi) . s {hi - a) 



a.2 e, f r cos cos sm sm 

Oi 020301 \ 04 ' 



71727303 / s(x-hi) 
OLi -TTTT ( COS sm sm cos 



(fe) 

(51^2^302 \ 04 

7i727303 / s{x- h^)^ '''^ 



0:2 cos sm cos Sill 

0^020301 \ a4 ' 

717273O4 / . S (x - /l3)y'^ 

012 f ^ ^, — sm sin cos cos 



> 

03 




S (/l3 - 


/l2) 


03 




,s (hs - 


/J2) 


03 




s(/l3 - 


/i2) 



O2 




S (/l2 - 


hi) 


02 




s(/l2 - 


hi) 


02 




S (/l2 - 


hi) 



Ol 




s {hi - 


a) 


Ol 




s {hi - 


0) 


Ol 




s(/ii - 


a) 


Ol 




s (/ii - 


0) 



^1(52(5303 \ 04 / as 02 Ol 

7I7273O204 / . s(a;-/i3)y'^ . s{h3-h2) . s {h2 - hi) . s {hi - a) 
a.2 „ I sin ^ sm sm sm ■ 



5^^2^30301 \ 04 / 03 O2 Ol 

717273O4 / . s{x-h-i)\^^^ s(/i3-/i2) . s(/(2-/'i) s{hi-a) 
— Oi2 f „ „ — sm cos sm cos ■ 



51^2^302 V / ^3 ^2 Ol 

7i727304 / . s{x-h-i)\ s (/13 - /12) s{h2-hi) . s {hi - a) 



- a2 -TTTTTi — sm COS cos sm 

"1O2O3O1 \ 04 / as 02 Ol 

+0 ( Isf'^ exp ['^1^203 {x — /13) + 010204 (/13 - h2) + 010304 (/i2 - hi) + 020304 {hi — a)] 



01020304 

(36) 

// a2 = 

(fc) oiai / . s{x-a)\^''^ ( k^2 \t\{x-a)\ 
4>i\ {x) = — ( sm — 1 +0 Ms| exp 1 , (37) 



,()=).. oiai 71 / s{x — hi)\^''^ . s{hi — a) 
(t>2\ {x) = — cos sm ■ ' 



s 5i \ 02/ Ol 

027I / . s{x-hi)\^^'' s{hi-a) 
— a2 — f7 sin — cos ■ 



oi^i V 02 / Ol 

+ O (\sf-' exp imai(^-M + a2(/^i-o)]'j ^^^^ 
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(fc) 717204/ s{x-h2)y' s(/i2-/ii) . s{hi-a) 
4>-^\ (x) = — Q!i — r-r- cos cos sin 



7172^2/ s{x-h2)Y ' . s{h2-hi) s {hi - a) 
ai — ^rr~ cos sm cos ■ 



(fe) 

sS[S2 y"" 03 / 02 ai 

71720103/ s (x - /i2)\ ''''^ . s{h2-hi) . s{hi-a) 
ai — —— sm sm sm ■ 



s5\5i2a2 \ CLa J 0-2 ai 

7i72«3 / . s(a;- /i2)V*'' s{h2-hi) s{hi-a) 
c«i — f , f , sm — cos— cos 



s5[5'2 \ 0,3 J 0-2 o-i 

{hi - a)] 

(39) 



+ O ( |s|''~^ exp 1^1 ["i"2 (a; - /t2) + aias (/12 - /ii) + 0203 (/ii - a)] 



010203 



(fc) 7172730-1/ s(a;-/i3)y'^ s(/i3-/i2) s(/i2-/ii) . s(/ii-a) 

(2^) = ~CKl e c c cos — COS — COS — SlU — 

Sdid203 \ O4 / O3 02 Oi 

71727302 / s{x-h3)\ ^''^ S (/l3 - ^2) . S (/l2 - hi) s {hi - o) 

— ai — COS COS ^ sm ^ cos ^ 

S0i020a, \ a4 J 03 02 oi 



7172730301/ s{x-h3)Y' . s{h3-h2) . s{h2-hi) . s {hi - a) 
+ ai — cos sm ^ sm ^ sm ■ 



(fe) 

S(5i^2^3a2 V O4 / 03 02 Oi 

7i727303 / s{x-h3)\^''^ . s{h3~h2) s{h2-hi) s {hi - a) 
— ai — ^rTrr~ cos sm cos cos 

SO1O2O3 \ 04 / 03 02 Oi 

7172730401/. s{x-h3)y'^ . s(/i3-/i2) s{h2-hi) . s{hi-a) 
+ ai — - , „ sm sm cos — sm ■ 



04 




s{x — 


M 


04 




s {x — 


h3) 



s5i525'^a3 \ ai J 03 02 Oi 

,.,-,0- - , • V s{h3-h2) . s{h2-hi) s{hi-a) 

+ ai — - sin sm sm cos 

sdi0203O3 \ 04 y O3 O2 Oi 

^'"^ s{li:>,-li2) . s{h2-hi) . s{hi-a) 
ai ' I sm — ^ — I cos sm sm ■ 



s6i6'26'3a2 \ tti J 03 02 oi 

717273O4 / . s{x-h3)\ ^''^ s (/13 - /12) s{h2-hi) s {hi - a) 
~ Q!i — „ sm cos — cos — cos ■ 



sd'i525'^ V 04 / 03 02 oi 

, ^ / I ifc-2 1*1 [010203 {x - /13) + 010204 {hs - /12) + 01O3O4 (/i2 - hi) + 020304 {hi - a)] 
+(J I \s\ exp 



01020304 

(40) 

Proof. We may prove these formulas using the similar method as [15, Lemma 
1.7 p. 9-10]. ■ 

Theorem 3.1. Let X = s"^ , t = Ims. Then the characteristic function co (A) 
has the following asymptotic formulas: 
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Case 1: If I3'2j^0, a2j^ 0, then 

M \ 3 a' 717273 . s{b- hs) s (/13 - /12) s (/12 - /ii) s {hi - a) 

ai0\020z a4 03 02 ai 

3o/ «27l7273 . 5(^-/13) s(/i3-/i2) . s {h2 - hi) . s {hi - o) 
— a2S P2 sill ^ cos ^ ^ sm — sin 



' 01045^^2^3 O-A fl3 O2 fll 

3o/ 03717273 . 5(6-/13) . s(/i3-/i2) . s (h2 - hi) s{hi-a) 
a2S P2 „ sm — sm — sm cos ^ 

0204010203 04 03 02 Ol 

3ol 03717273 . 5(6-/13) . s(/l3-/l2) S (/12 - /ll) . S (/li - o) 

- Q!2S P2 „ „ c sm ^ sm ^ cos ^ ^ sin ^ 

01040^1^0203 O4 O3 02 Ol 

3o/ 717273 5(6-/13) . s{hs-h2) .s(/(2-/'i) .s(/ii-o) 
+ a2S P2 — „ cos sm cos cos 

O3O1O2O3 04 03 O2 Ol 

3n, 027^7273 5(6-/13) . 5(/i3-/i2) . s{h2-hi) . s{hi-a) 
+ a2S^P2 „ e- „ cos sin ^ sm ^ sin ^ - 

OiO30!|^0203 04 03 02 Ol 

3 0/ 717273 S (6 - hs) 5 (/13 - /I2) . 5 (/l2 - /ll) 5 (/ll - o) 

+ a2S P2 TTTTT cos COS Sm COS 

02C)1C)2C'3 04 O3 O2 Ol 

3 nr I1I2I3 s{b- /13) 5 (/I3 - /I2) 5 (/l2 - /ll) . 5 (/ll - o) 

+ a2S P2 — f, f, „ COS ^ COS ^ ^ COS sm ^ 

aid[d2d^ 04 03 02 Ol 

+0 ( |5|^ exp 1*1 [<^i02a3 (6 - /13) + 010204 (/13 - /12) + 0103O4 (/i2 - hi) + 020304 (/ii - o)] 

01020304 

(41) 

Case 2: If (5'2T^0, a2= 0, then 

IW 2/j/Ol7l7273 . 5(6-/13) 5(/l3-/l2) 5(/l2-/ll) . 5 (/ll - o) 

a;4 (A) = — ais P2 — sm — ^cos— ^cos— ^sm- 



04^152^3 CI4 03 O2 ffll 

2a/"27l7273 . 5(6-/13) s(/l3-/l2) . S{h2-hi) 5 (/ll - o) 

ais P2 — sm ^ COS ^ sm ^ cos 

04010203 04 03 02 Ol 



2 „, 0103717273 . 5 (6 - /I3) . 5 (/I3 - /I2) . 5 (/I2 - /ll) . 5 (/ll - o) 

— ais P2 -7, sin — sm — sm ^ sin — 

020401021)3 O4 03 O2 Ol 

2 0/ a37l7273 . 5 (6 - /13) . 5 (/I3 - /12) 5 (/l2 - /ll) 5 (/ll - o) 

— ais P2 — sm — sin ^ cos ^ cos ^ 

040^0203 04 03 02 Ol 

2 0/0171727.3 5(6-/13) . s(/l3-/l2) 5(/l2-/ll) . 5 (/li - o) 

— ai5 P2 — ^ ^ ^, COS sm cos sm 

O3O1O2O3 04 O3 ©2 Ol 

2 Ql «27l7273 S{b-h3) . S (/13 - /12) . 5 (/l2 - /ll) 5 (/ll - o) 

+ ais P2 — COS ^ sin — sin ^ cos — 

O3O1O2O3 04 03 O2 Ol 

20/ "^1717273 s{b-hs) 5(/i3-/i2) . s(/i2-/ii) . s{hi-a) 

— ais P2 — cos cos sm sm 

O2O1O2O3 04 O3 O2 Ol 

20/ 7i 72 73 5(6-/13) 5(/i3-/i2) s{h2-hi) s{hi-a) 
+ ais P2 „ „ „ cos — cos — cos cos — 

0^^0203 O4 03 02 Ol 
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+o(\s\ex 1^1 [Q'i^2Q3 {b - hs) + 010204 {hs - /12) + 0103^4 {h2 - hi) + a2a3a4 {hi - a) 
\ 01020304 

(42) 

Case 3: If P'2=0, a2 0, then 

rw 20/ 717273 5(6-/13) s(/i3-/i2) s{h2-hi) s{hi-a) 

uji (A) = a2S Pi ^ ^ ^ cos cos cos cos 

010203 04 O3 O2 Oi 

2o/'^27i7273 5(6-/13) s{h3-h2) . s{h2-hi) . s{hi-a) 

— a2S Pi — - - cos ^ cos ^ sm — sm 

aio[o20s 04 03 02 oi 

2^/03717273 5(6-/13) . s(/i3-/i2) . s{h2-hi) s{hi-a) 

— a2S Pi — ^ „ ^ cos sin — sin — cos — 

a20ic>2<^3 a4 da 0-2 ai 

2 pi a37i7273 s{b-h3) . s (/13 - /12) s (/12 - /ii) . s {hi - a) 

— a2S Pi — - cos sm cos sm 

0105^02^3 O4 O3 O2 Oi 

20/^4717273 . s{b-hs) . s(/i3-/i2) s{h2-hi) s{hi-a) 

— a28 Pi — r^TT ^ sm ^ cos ^ cos ^ 

03010203 04 03 02 Oi 

2 01 02047^7273 . S{b-h3) . 5 (/13 - /12) . 5 (/l2 - hi) . 5 (/li - O) 

+ a2S Pi Frrrr sm sm sm sm 

030105^0203 04 03 O2 Oi 

2 „, 047x7273 . 5(6-/13) s(/l3-/l2) . s(/i2-/ii) s {hi - a) 

— a2S Pi — TTTTT" sm ^ cos ^ sm — cos ^ 

O2O1O2O3 04 03 O2 Oi 

2 0/ 047l7273 . 5 (6 - /13) 5 (/13 - /12) s{h2- hi) . 5 (/li - o) 

+ q;2S Pi — „ „ J sm cos cos sin 

O1OXO2O3 04 03 02 Oi 

+0 (\s\ exp 1^1 [010203 (6 - /13) + 01O2O4 {hs - /12) + 01O3O4 (/i2 - hi) + 02O3O4 {hi - 
\ 01020304 

(43) 

Cose ^ : 7/ /32 = 0, a2 = 0, then 

„/ 01717273 5(6-/13) s(/i3-/i2) s{h2-hi) . s(/ii-o) 

^4 (A) = — Qi5Pi cos COS COS sm 

010203 04 03 02 Ol 

a, 027(7273 s{b- /13) s (/13 - /12) . 5 (/l2 - /li) 5 {hi - a) 

— aiSPi — :7-r-j COS — COS — Sm COS — 

O1O2O3 O4 O3 O2 Ol 

0103717273 5(6-/13) . 5(/l3-/l2) . S {h2 - hi) . 5 (/li - o) 

— aispi — COS sm sin sm 

O2O1O2O3 04 03 02 Ol 

a/037i7273 5(6-/13) . s{hs-h2) s(/i2-/ii) s{hi-a) 

— aispi ^, ^ — cos sin cos cos 

0;^02O3 04 O3 02 Ol 

-,0104717273 . 5(6-/13) . 5(/l3-/l2) S{h2-hi) ^ "'^ 

— aispi - sin sin cos sm 

03O1O2O3 04 ~ 



o 



f "-z""! 1 z n ■ " \" '"a) . " ya '"z) . " yz '"l/ 

■ aispi sm sm sm cos 

O3<Ji<J203 04- 



03 




s{hs - 


h2) 


03 




5(/l3 - 


h2) 



02 




5 (/l2 - 


hi) 


O2 




5 (/l2 - 


hi) 



Ol 




5 {hi - 


a) 


Ol 




5 {hi - 


a) 



I 1*11^11 ,1 ,2 (3 . " V" '"'3/' " yo '"z/ . " yz '"Lj . 

— aispi sin cos sin sin 

O2did203 O4 O3 O2 Ol 

o/047i7273 • 5(6-/13) s{hs-h2) s{h2-hi) s{hi-a) 
+ aispi sm cos cos — cos 

0^0203 04 03 02 Ol 
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/ ^ 1^1 [aia2a3 (6 - /13) + 010204 (fes - ^2) + aia3a4 (/t2 - hi) + 0203^4 (/ti - a)] ' 
V 01020304 

(44) 

Proof. The proof is completed by substituting (36) and (40) into the represen- 
tation 

W4 (A) = A [I3[(f>ix (b) - (i'^cj>'^^ (6)] + [/3i04A [b) - M'^^ (&)] 

= -A/3^0;A {b) + A/3;04A (6) + /5l04A (6) - /?2</'4A (^) • (45) 

■ 

Corollary 3.1. The eigenvalues of the problem (l)-(9) are hounded below. 
Proof. Putting s = it {t > Q) in the above formulas, it follows that a;4 (— t^) — >■ 
00 as t ^ 00. Therefore, oj^ (A) ^ for A negative and sufficiently large. ■ 

4 Asymptotic formulas for eigenvalues and eigen- 
functions 

Now we can obtain the asymptotic approximation formulas for the eigenvalues 
of the considered problem (l)-(9). 

Since the eigenvalues coincide with the zeros of the entire function a;4 (A), 
it follows that they have no finite limit. Moreover, we know from Corollaries 
2.1 and 3.1 that all eigenvalues are real and bounded below. Hence, we may 
renumber them as Ao < Ai < A2 < listed according to their multiplicity. 

Theorem 4.1. The eigenvalues A„ = s^, n = 0,1,2,... of the problem 
(l)-(9) have the following asymptotic formulas for n — >■ 00 : 

Case 1: /f 7^ 0, a2 ^ 0, then 

0ia2a3047r in — 1) 

" 010203 (6 — hi) + 010204 (/13 — h2) + 01O3O4 (/i2 — hx) + 020304 (/ii — o) 
+ 0(i). (46) 

Case 2: If ^'2 7^ 0,a2 = 0, then 

0ia203047r {n — ^) 

" 01O2O3 (6 — /13) + 01O2O4 (/i3 - /12) + 01O3O4 (/i2 — hi) + 020304 {hi - a) 
+ 0(i). (47) 

Case 3: // = 0, a2 ^ 0, then 

010203047: (n — ^) 

" 01O2O3 {b - /13) + 01O2O4 (/13 - /12) + 01O3O4 (/i2 - /ii) + 02O3O4 {hi — a) 
+ 0(i). (48) 
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Case 4- = 0, a2 = 0, then 



aia2a3a4TTn 



010203 (6 — hs) + 010204 {hs — h2) + 010304 (/i2 — hi) + 02O3O4 [hi - a) 

(49) 



.oil 

Proof. We shall only consider the first case. The other cases may be considered 
similarly. Denoting ojq (s) and UI2 (s) the first and O-term of the right of (42) 
repectively, we shall apply the well-known Rouche's theorem, which asserts that 
if / (s) and g (s) are analytic inside and on a closed contour C, and |.g(s)| < 
1/ (s)| on C, then / (s) and f (s) + g (s) have the same number zeros inside C, 
provided that each zero is counted according to their multiplicity. It is readily 
shown that |ZUi (s)| > \iJ2 {s)\ on the contours 



Cn — <S G C||s| 



7raia20304 (n + i) 



O1O2O3 {b - /13) + 01O2O4 (/i3 - h2) + 010304 (/i2 - hi) + 020304 {hi — a) 



for sufficiently large n. 

Let Ao < Ai < A2 < ... be zeros of w (A) and A„ — s^. Since inside the 
contour C„, ZJi (s) has zeros at points s = with multiplicity 4 and 

fc7raia2a304 



s = 



2 [010203 {b - /13) + 01O2O4 (/13 - /12) + 010304 (/i2 - /ii) + 02O3O4 {hi - a)] ' 
fc = ±l,±2,...,±n 



with multiplicity 1, and the number of zeros is 2n + 4, it follows that 

7701020304 (n — 1) 

s — ^ 

" 010203 {b - /13) + 01O2O4 {ha - /12) + 01O3O4 (/i2 - hi) + 02O3O4 {hi — a) 

where 6n = 0{1), more precisely 

1^ I ^ 7701020304 

" 4 [010203 (6 - /is) + 01O2O4 (/13 - /12) + 010304 (/i2 - hi) + 02O3O4 (ft-i — a)] 

for sufficiently large n. By substituting this in (41), we derive that (5„ = O (i), 
which completes the proof. ■ 

Recalling that (f)(x,\n) an eigenfunction according to eigenvalue A„, by 
putting (46)- (49) in the (33)- (40) for fc = 0, we derive the eigenfunctions: 



In case 1 



'(a;, A„) = < 



a2 cos 



a2a3a4(x— a)7r(n— 1) 



Si 



■ COS 



010203(6— /i3)-|-ai 0204 (/i3 — /»2)-|-aia3a4(/(2 — /(I )-|-a203a4(/ii— a) 

+0(i), xe[a,hi) 

[ai 0304 (2: — /ii)-t-a2 03a4 (/ii — a)]7r(n— 1) 



i5ii52 



COS 



a27i7273 (.Qg 
<5i52<53 



010203(6— ^3) +010204 (?i3 — /t2)-l-ai£i3 04 ('12 — 'Jl)-|-a2a304(ftl — o) 

+ 0(i), X€{hi,h2) 

[aia2 0,3{x — h2)-\-aio,2a4{h2 — hi)-{-a2a3a4{hi—a)]7T{n—l) 
aia2a3{b—h3)-\-aia2a4{h3 — h2)-\-aia3a4ih2 — hi)-\-a2a3a4ihi—a) 

+0(i), XG{h2,h3) 

[ai0204(x — /t3)-F0l 0203(^3— ^2)-|-Ol0304(fc2—fel)-|-02 0304(/n—o)]7r(7t—l) 
010203(6 — ft3)-|-0i0204(/»3—/l2)-F0l0304(/t2 — /n)-|-02a3a4('H—a) 

+0(i), x€{hs,b] 
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In case 2 



{x, Xn) 



In case 3 



(j){x, Xn) 



X sm 



^ 010203(6— '»3)+aia2a4('J3 — 'i2)+ai 0304 (/i2 — /ii)+a2a304('»i— a) 

—cxi 7 — ri 



020304 TT^^ 

020304 (a: — a) TT^^n — ^ ^ 



+ 



X sm 



0102 03(fe — /l3) + 01 02 04(^13— /l2) + 01 03 04(712— /ll) +02 0304 (/ll—o) 

X e [a, hi) 

ai7l 010203(6— fc3) + Ol0204(fe3 — fa2)+Ol 0304(^2— fcl)+020304(/n—o) 

^1 0203a47r(n— 

[010304(2; — /ii)+a2 03 04(/ii— o)]7r^n— 



+ 



X sm 



+ 



X sm 



010203(6- /is) + 01 0204 (/l3—/i2) + 01 0304 (/l2—/ll) +02 0304 (/ll—o) 

X e {hi, /12) 

ai7l72 010203(6 — fe3)+Oi 0204(^3 — ft2)+Ol03 04(/t2— fcl)+0203 "4(^1— o) 
'5l^2 02O3O47r(n— 5) 

[010203(2; — /l2) +01 03 04(/l2—^l) + 02 03 04(/li—o)]7r^n— 
010203(6- /is) + 01 0204 (/i3—/l2) + 01 0304(712-/11 ) +02 0304 (/j-l-o) 

X G (/l2, /is) 

ai7l7273 Oi0203(6-/t3)+Oi0204(/t3— /l2) + Oi03 04(/(2— /ll)+0203a4(/li-o) 
l5l52l53 2 3 47r(n— ^) 

010204(2; — /13) +01 02 03(^13- /12) +01 0304(712- /ii) +02 03 04 (/i 1—0)] TT^n—^^ 

0102 03(/> — /is) + 01 0204 (/13— 712) + 01 03 04(712— Til) +02 03 04(6-1-0) 

X e {hs, b] 

02 0304(x — o)-7r(n— ^) 



a2 COS 



Si 



■ COS 



Hmcos 

"102 



0:2717273 (.Qg 

<5l5253 



010203(6— /l3)+Ol 0204(713— 62)+aia304(/l2 — Til )+a2 0304(/n—o) 

+0(i), xe[a,/ii) 

[aia3a4(a; — /ii) +0-2 as 04(^1 — a)]7r(^n— 
tiiti2a3(^'— ^3)+aia2a4(/i'3 — 'i2)+aia3a4(/j'2— /n)+a2a3a4(^i — o) 

+ 0(i), XG(/li,/l2) 

[0l0203(x — 7l2) + Ol 0304 (/l2— Til )+02 0s04(/ll—o)]7r(^n—ij 

010203(6— 7i3)+oi 0204(713 — 7i2)+aia304(/i2 — Til )+a2 0304(711— o) 

+ 0(i), Xe(/l2,/l3) 
[aia2a4(x — ^;-3)+a'ia2a3(/i3 — ^2)+aia3a4(/i2 — /ii)+a2a3a4(^i — a)]7r(n— "I) 
010203(6— /i3)+aia2a4(/i3 — ^2)+aia3a4(/j'2 — /ii)+a2a3a4(/ii— a) 

+0(i), a;e(/i3,6] 
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In case 4 



(j) (x, A„) 



-ai 



aia2a3(b— fc3)+aia2 04(^3 — ^2)+aia3a4(fe2 — /ti)+a2a3a4(fci — a) 



X Sin 



X sm 



X sin 



020304 Trn 
020304(3; — a)7rn 



010203(6 — ?l3)+Ol02a4(/l3—/l2)+Oi0304(?l2 — '»l)+02a304(?ll— a) 

.r G [a. hi) 

ai7i ai».2H;i(/>-/';i)— »ia2" l(/';i-/'2)+"a";l"a(/'2-/'l)+"2";i" l(/'l -o) 



Si 



[aiaza4^{x — hi)-\-a2a'ZaA{h\—a)\-nn 



0\ 



010203(6 — ?13)+0102 04(/t3 — ?i.2)+Oia3a4(?l2 — 'tl)+a20304(ft.l— a) 

X e (/ii, /12) 

ai7l72 010203(6— ^3) +0102 04 (^3 — ^2) +01 0304(^2— ^l)+020304 (hi— g) 
(5i(52 a203a4'7rn 



[010203(3: — /t2)+oi 0304(^2 — hi )+0203 04(^11- a)]7rn 



X sm 



010203(6 — h3)+Ol02 04(/l3 — /l2)+Oi03a4(?l2 — hl)+a20304(/M— o) 

ai7i7273 010203(6— /i3)+ai0204(/t3 -h2)+oi03 04(/t2 -hi )+02a304 (^1-0) 

S1S2S3 02 03 04 7rn 



[010204(3: — hs) +01 02 03(^3 — ^2) +01 0304(7^2 — hi) +02 03 04 (h 1—0)] 7rn 
010203(6 — h3)+Ol 02 04(h3—h2)+Ol 0304(h2 — hi )+0203 04 (hi— o) 

a; G (/i3,&] 



All these asymptotic formulas hold uniformly for x € [a, /ii)U(/ii, /i2)U(ft.2, /i3)U 
ih3,b]. 
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